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Many models have been studied that contain more than one species of dark matter and some 
of these couple the Cold Dark Matter (CDM) to a light scalar field, fn doing this we introduce 
additional long range forces, which in turn can significantly affect our estimates of cosmological 
parameters if not properly accounted for. ft is, therefore, important to study these models and their 
resulting cosmological implications. We present a model in which a fraction of the total cold dark 
matter density is coupled to a scalar field. We study the background and perturbation evolution and 
calculate the resulting Cosmic Microwave Background anisotropy spectra. The greater the fraction 
of dark matter coupled to the scalar field and the stronger the coupling strength, the greater the 
deviation of the background evolution from ACDM. Previous work, with a single coupled dark 
matter species, has found an upper limit on the coupling strength of order O(0. f). We find that 
with a coupling of this magnitude more than half the dark matter can be coupled to a scalar field 
5-H ' without producing any significant deviations from ACDM. 
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I. INTRODUCTION 

o. 

Since the first discovery that visible matter does not account for all the mass in the Universe [1[ astronomers 
have known that there is something missing from our basic understanding of the Universe. The two most popular- 
explanations for this missing mass have been that there is some other matter present that we cannot see, in the 
form of Cold Dark Matter (CDM), or that our theory of gravity is inadequate, giving rise to a plethora of modified 
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theories of gravity (for a review see e.g. Ref. |2|). More recently, many more accurate measurements of the cosmological 
parameters have been made e.g. Refs. [J-|g|, constraining the amount of CDM in the universe to be roughly 22% of the 
total energy budget. Cosmological data analyses typically assume a ACDM model in which the accelerated expansion 
is due to a small cosmological constant and gravity is described by General Relativity, universally coupled to all species 
of matter (visible and dark) . Any deviation from these assumptions could affect the best fit cosmological parameters 
(see e.g. Ref. p}). In this paper, we will alter the latter assumption, allowing for violations of the Equivalence Principle 
in the dark sector. In particular, we will study a generalised Brans-Dicke scenario in which some species of dark matter 
couple to the Brans-Dicke scalar, while others do not. ACDM can be identified with the limit in which the coupling 
vanishes, or else the limit in which the fraction of dark matter coupled to the Brans-Dicke field goes to zero. Here 

p^i | we will explore what happens as the coupling and the fraction of coupled dark matter are increased. 

Since the exact nature of dark matter is unknown, it is important to study different models and their effects 

P^L upon cosmological observables. In particular, plenty of recent work has focussed on the case where dark matter is 

coupled to a quintessence scalar field playing the role of dark energy e.g. Refs. [a-Uo|. Because the quintessence 

^-H ' field is necessarily light, this will mediate an additional long range force on dark matter, breaking the gravitational 
equivalence with baryonic matter. The astrophysical and cosmological implications of an additional long range force 
acting exclusively on dark matter were studied long ago by Gradwohl and Frieman |Tl| - fl3l |. and later by Gubser and 
Peebles [1J, [l5|], and others [P6h18| . Prior to that, effects on time varying Newton's constant were used to place 
bounds on a generalised Brans-Dicke scenario in which visible and dark matter both couple to the Brans-Dicke field, 
but with differing strengths [19j . 

If one assumes more than one species of dark matter, as suggested by the richness of string compactifications [2fJ, [21| , 
one can also entertain the possibility that the Equivalence Principle might be violated within the dark matter sector 
itself. Indeed, there is no good theoretical reason to maintain it once we have abandoned the universal coupling 
between matter and gravity by allowing dark matter to feel an additional long range force. This possibility was 
explored to some extent in Ref. [22], where different dark matter species are coupled to a scalar field with different 
strengths. In that work the scalar field acts under the influence of a quintessence potential, playing the role of dark 
energy. This gives rise to chameleon-like behaviour [23|, |2J] , with the scalar settling in to a minimum of the effective 
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potential controlled by the density of the coupled fields. In contrast, here we assume that dark energy is a cosmological 
constant, and that there is no bare potential for the scalar - in other words, it is exactly massless. This alters the 
picture considerably as there is no longer a local environmental minimum for the effective potential. Such a scenario is 
easily achieved in higher dimensional compactihcations. For example, in the two-brane Randall-Sundrum model 25], 
the warped extra dimension gives rise to exactly the sort of scenario we are interested in, with an exactly massless 
scalar (the radion) coupled with different strength to matter on different branes [2a ]. 

If they are not properly accounted for, additional long range forces in the dark sector can dramatically affect our 
estimates for cosmological parameters. This point was emphasized in Ref. [27| . where ignorance of the additional long 
range force was shown to result in an erroneously inferred equation of state for dark energy 1 . We can illustrate the 
principle by considering the orbit of a planet around a star. One could compute the mass of the star by measuring 
the acceleration of the planet, and then use Newton's law of gravitation to infer the mass. However, if there were an 
additional long range force acting between the star and the planet, and we were ignorant of this, it is clear that we 
would infer the wrong value for the star's mass. For cosmology, the lesson we learn is that we should ask whether or 
not the relevant observational probe is sensitive to the additional long range force. When this only acts on (a fraction 
of) dark matter, the supernova probe will be directly sensitive to the force as supernovae reside in galaxies, which are 
dominated by dark matter. Galaxies therefore follow accelerated rather than geodesic trajectories and this ought to 
be properly accounted for when analysing supernova data [271 ] - In this paper, we will focus on a more indirect probe 
of the scalar force - the Cosmic Microwave Background (CMB). CMB photons are part of the visible sector, so they 
do not feel the extra force directly. However, the CMB temperature fluctuations are sensitive to density fluctuations, 
which are affected by any additional long range forces acting on any form of matter, be it visible or dark. 

The paper is structured as follows: in Sec. [H]we describe the setup for our model, then in Sec. [HI] we study the 
background evolution and examine the extent to which it deviates from ACDM. In Sec. IIVI we study the evolution 
of perturbations in our system and the effect the coupling has on CMB anisotropies. In addition to comparing to 
standard ACDM, we discuss the physics behind the deviations and find that the main physical effects come from the 
increased density of CDM at early times, a result that is only indirectly related to the coupling. Finally, we conclude 
in Sec. El 

II. SETUP 

Let us begin by describing our setup. We consider a generalised Brans-Dicke theory with a cosmological constant 
responsible for late time acceleration of the universe. The matter content is divided into two components, one of which 
consists of baryons, radiation and some fraction of dark matter, and the other which is made up of the remaining 
fraction of dark matter coupled to the massless scalar field. The action can be written in the Einstein frame as follows 
(see also Ref. [H) 



S = d 4 x 
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+ S SM [g^, ■•■] + S c [g^, ...] + S q \p^,ip q ] , (1) 



where g is the determinant of the metric g llVl R is the corresponding Ricci scalar, V is the corresponding covariant 
derivative and (f> is the scalar field mediating the additional long range force. The action SsmIs^i/j ■••] contains baryons, 
photons and massless neutrinos and Sclg^, ...) contains uncoupled CDM. The ellipses denote the standard model fields 
in S'sm and uncoupled dark matter field in S c . None of these fields couple directly to the scalar. The action Sqlg^, ip q ] 
describes dark matter fields that do couple directly to the scalar. Indeed, coupled dark matter particles, Sifjq, follow 
geodesies of the conformally related metric, 

g„u = e 2Q V„ , (2) 

where a is the dimensionful coupling constant, which, without loss of generality, we assume to be positive. The field 
equations describing the system are given by 

G» + AS; = SvrG (T£ M> + If c)v + Tfc v + % v ) , (3) 



1 In Ref. 12711 the additional long range force is mediated by a gauge field and yields an additional repulsive force. Here we are focussing 
on a scalar mediator which will give rise to an additional attractive force, but it is clear that the principles discussed in Ref. [27l | will 
still apply, although the effects will work in the opposite direction. 



where □ = ^"VuV,, is the covariant d'Alembertian, 
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are the energy-momentum tensors for standard model fields and uncoupled CDM, and 

r ( ; )M = v>v^-^(v^) 2 , (6) 

is the energy-momentum tensor for the scalar. In addition, we define 

*?,)* = «***K>„ ' ( 7 ) 

where the energy-momentum tensor for the coupled fields is covariantly conserved in the conformally related Jordan 
frame, and is given by 



T ^ = ~7^/ a s^- (8) 

Note that T/\ is not covariantly conserved in the Einstein frame, but satisfies V M T/\ = aT?-. V„(/>. 

III. BACKGROUND COSMOLOGY 

To leading order, we consider a homogeneous and isotropic universe described by a flat, Friedmann-Robertson- 
Walker metric, ds 2 = — dt 2 + a 2 (i)dx 2 , and a homogeneous scalar <f> = <fi(t). Each matter field contributes a perfect 
fluid 

T (sm)i/ = di&g(-psM(t),psM(t),PsM(t),psM(t)), Tfa v = diag(-p c (t), 0, 0, 0) 
T^ = diag(-p,(«), 0,0,0) (9) 

where psm = Pb + Pi and psu = Pb + Pr, where b denotes the baryons and r collectively denotes radiation and 
neutrinos. We have that pb = 0, p r — p r /3 and p^ — p^ — 4> 2 /2. Now, neither Tl\ nor T/\ is covariantly conserved 

in the Einstein frame, and by association, neither p q nor p q = e ic "^p q satisfy the usual conservation law for matter 

(i.e. d(a 3 p q ) / dt ^ 0, d(a 3 p q ) / dt ^ 0). However, we can use the fact that T/\ is covariantly conserved in the Jordan 

frame to show that p* — e 3 °"^ p q = e~ a ^p q scales like 1/a 3 [23|, [24j. The background evolution equations are given by 
the Friedmann equation, 

H 2 = ^ (p A + psm + p c + \<P 2 + p*e a ^j , (10) 

and the homogeneous equation of motion for the scalar 

$ + 3H$ + ap*e a * = 0, (11) 

where the Hubble parameter H = a/a and overdots represent differentiation with respect to cosmic time t. Note that 
Pa = A/87rG is constant, while p\, oc p c oc p* ex 1/a 3 , and p T oc 1/a 4 . It is also convenient to define the following 
density fractions 

_ 87rGp A _ 8nGpn „ _ 4ttG0 2 _ 87rGe°"V. _ 87rGp r 

"A- 3R2 , i«„- 3ij2 , »*- 3i/2 , »,- 3ff2 , Sir- 3jff2 

where n denotes either baryons or uncoupled CDM, and in analogy with quintessence models, we have identified 
£!de = ^A + Q(j>- We note that this grouping of dark components is somewhat arbitrary [28j |. 

Let us now consider the resulting evolution of the cosmological background. The coupling to (some) dark matter 
endows the scalar with an effective potential V e fi(<l>) = p*e Q *. 

Unlike the standard chameleon scenarios |23l |24| where a local minimum exists at finite <f> (due to the existence of 
a bare potential for 0), our effective potential is of run-away form, with a minimum only at +oo. The only way the 
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FIG. 1: The evolution of the density parameters in our model (thick lines) compared with ACDM (thin lines), for varying 
coupling strengths, a, and varying proportions of coupled dark matter (CCDM) such that CCDM=4% means that 4% of the 
total energy density in the universe is made up of coupled dark matter. The energy density of coupled dark matter is £l q , 
fl c is the energy density of uncoupled dark matter and the total dark matter energy density is fic to t — Q q + Q c . The energy 
densities have been normalised to their present day values i.e. the system was evolved to give the same energy densities today 
as in ACDM. The coupling constant a is in units of 1/M P \, where M p \ is the Planck mass. 



scalar can be held approximately fixed is by Hubble friction, and only when 2 H ^$> a 2 p*. Therefore it is difficult to 
imagine that the scalar field does not significantly affect the background evolution (unless a <C 1/M p i) in contrast to 
the assumptions made in Ref. [12| where the background is always taken to be equivalent to ACDM. 

These expectations are born out by our numerical solutions plotted in Fig. [T] Here we show the background 
evolution for various fractions of coupled dark matter today and various coupling strengths, a. We compare these 
to the corresponding evolution in a ACDM universe, assuming ACDM best-fit cosmological parameter values. We 
normalise each coupled dark matter model such that they share the same baryon, total dark matter, and dark energy 
density fraction today (Q£, f2° , Q^e) as the best fit ACDM model. In the limit a — 0, we recover ACDM where 
rip E = tt\ and Q° = SlJ?. This normalisation is achieved by performing a binary search over the initial condition fa 
and over the cosmological constant density p\ to find the correct initial conditions necessary in order to realise the 
desired background. We set fa — since we begin the evolution deep inside the radiation dominated epoch, where we 
should expect the field to be frozen by Hubble friction. If the coupling is large enough, we can easily see the effect of 
the evolving scalar field as we move back in time. Initially, at low redshift, the main effect is to decrease the density 
fraction stored in baryons (and dark matter), relative to ACDM. This is to be expected because the mutual coupling 
to gravity allows energy to pass from matter to the kinetic energy of the evolving scalar. At higher redshifts, Hubble 
friction suppresses the kinetic energy of the scalar, but we still see the effect of the scalar's previous evolution on 
the relative density of radiation to matter. In particular, matter-radiation equality happens at higher redshift in our 
coupled dark matter scenario. This is because we need more matter at high redshift to compensate for the relative dip 



2 This can be seen by dropping the acceleration term in Eq. <1H , and rearranging to give a 2 p t /3H Ri d(e °"t')/dt. This derivative is then 



in the matter density fraction at lower redshift, if all cosmological parameters are normalised by their values today. 
This effect also manifests itself through a deficit in the amount of radiation at high redshift, compared with ACDM. 
In any event, it is clear that the background evolution is sensitive to coupling (some) matter to a light scalar. This 
effect is often neglected in the literature (see e.g. Ref. [Da. l29l. |3G| ) . 

To further analyse the background dynamics it is instructive to transform our system into autonomous form and 
investigate the presence of any late time attractor solutions. We use the standard definitions j31| along with an extra 
term to include the conformal factor, 
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The flatness condition gives the constraint 



Q b+C = I - x 2 - y 2 - z 2 - s 2 , 



(12) 



(13) 



where we have grouped baryons and the uncoupled CDM together as they both scale like a~ 3 . In these equations p 1 
and pa are the radiation and cosmological constant energy densities respectively. Differentiating these equations with 
respect to N = In a, which we denote using ', yields the following set of equations 

, (W \ 3 2 
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where 



§■ = - l -(3x 2 +y 2 -3 S 2 + 3). 
The system fEqs.[T4"|) generically admits 6 critical points at real values, listed in Table fl] 
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TABLE I: Critical points for the autonomous system (Eqs.[l4|. 



Three of these fixed points are familiar from ACDM: uncoupled matter domination, radiation domination, and A 
domination. As with ACDM our stability analysis reveals that A domination is an attractor, whilst the other two 
are saddles. This is to be expected from the cosmic no hair theorem [32j |. In our generalised Brans-Dicke model, 
we generically have three additional fixed points. The first of these is always present, and corresponds to an era in 
which the kinetic energy of the scalar dominates the cosmic evolution. This fixed point is unstable. Whether or not 
the system evolves to this fixed point depends on the choice of initial conditions. For |a| < y3/2 there exists a 
saddle point corresponding to an era in which the coupled DM component and the kinetic energy of the scalar scale 
with one another, dominating the expansion. For |a| > \/l/2 there is another saddle point, corresponding to a three 
way scaling between coupled DM, radiation, and the kinetic energy of the scalar. These properties are explicitly 
demonstrated in the phase-space plot shown in Fig. [51 where the evolution inevitably points towards the unique late 
time attractor corresponding to A domination. 



Phase Portrait in s-x plane showing A Domination 




FIG. 2: Phase portrait in s — x plane which shows that all trajectories eventually end in A domination. The red points are 
equilibrium points and the green are unstable fixed points. 



IV. PERTURBATION EVOLUTION 

Now that we have studied the background evolution of our system, we can move on to the evolution of perturbations 
and hence examine how the coupling affects structure formation and the CMB anisotropy spectrum. We will work in 
the synchronous gauge and to first order (see e.g. Refs. [3J, J34|) where the line element is 



As 2 



a(r){-dr 2 + (% + h i j)dx i dx j } , 



(16) 



and the metric perturbation hij can be decomposed into a scalar, vector and tensor part. Note that we have switched 
to conformal time r for convenience. Focussing on scalar perturbations, in Fourier space k we set 



K, = 1m, 



* 3 ~o*i>3 



(i, 



(17) 



and can use the following definitions 
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where each species n has coordinate velocity tv n )j, density contrast 5 n and velocity field gradient 6 n . We write the 
perturbed energy-momentum tensor for uncoupled CDM as 



6lf c)0 = -Sp c , 



«?-)« 



PcV(c)i 



-ST, 



(c)0, 



ST, 
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(19) 



where the uncoupled CDM pressure and anisotropic stress have been set to zero. Conservation of the CDM energy- 
momentum tensor implies 



SL 



■iv. 



(20) 



where here, and subsequently, primes denote differentiation with respect to conformal time. As is well known, there 
is a residual gauge freedom which allows us to consistently set 9 C = 0. The perturbations of the Standard Model 
components are treated in the usual way by evolving the full Boltzmann hierarchy in CAMB. 
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FIG. 3: CMB anisotropy spectra for various coupling strengths and proportions of coupled cold dark matter, CCDM. 



The coupled dark matter species undergo a non-standard evolution. Assuming that the perturbed contribution to 
the energy-momentum tensor in Einstein frame is given by 



5T ( q )0 = Sp q , 8T£ q)i = p q v (q)l = -5T l {q)0 , 



we can use V U T/^,, = aTj 1 -. V v <f> to demonstrate that 



«(,)i=0, 



(21) 
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9 q - -ti + a5<p' , 



6' = -9qH + a{k 2 54> - 4>'6 q ) , 



(22) 
(23) 



where S<p is the perturbation in the scalar field and % — a' /a. 

Finally, the evolution of perturbations in the scalar field is governed by the perturbed scalar equation of motion, 
and is given by 



6cf)" + 2H6<p' + k 2 6<p + -h'(f>' = -aa 2 5p q . 



(24) 



We now modify CAMB [351 ] to incorporate perturbations in these new species along with the background equations, 
Eqs. (ITU1) and (ITT1) . to calculate the CMB anisotropy spectra. As discussed in Sec. IIII1 we must evolve our back- 
ground first, before evolving the perturbations. To do this, an adapted version of the quintessence module available 
for CAMB [36| was used since this was designed to evolve quintessence models, which also do not have a ACDM 
background. 

In Fig. |3] we plot the resulting CMB anisotropy power spectra for various proportions of charged dark matter and 
coupling strengths, along with the spectrum for ACDM. The most interesting thing to note from this plot is the fact 
that even at a level of 16% coupled CDM with coupling strength a — 0.1, we still obtain a CMB power spectrum that 
is extremely close to ACDM. As both the coupling strength and proportion of CDM coupled increase, we see greater 
deviation from ACDM, to the point that a coupling strength of a — 0.5 or above is unlikely. These deviations are 
characterised by the large increase in power due to the Integrated Sachs- Wolfe (ISW) effect at low multipoles, and the 
lowering of the CMB peaks. The former has also been seen in similar models |37l.[38| and is mainly due to the late time 
variation of the gravitational potentials, the wells growing deeper in time as the scalar evolves to larger values. The 
latter is due to the fact that there is an increase in the coupled CDM density at early times, as is evident from Fig. [1] 
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This results in fewer baryons at the time of recombination, which in turn raises the sound speed and suppresses the 
amplitude of oscillations. Furthermore, if matter-radiation equality is earlier then the potential at the peak is deeper 
due to less pressure damping resulting in a large redshift as the photons climb out after last scattering, again, implying 
a smaller final anisotropy (for an excellent discussion of these effects, see Ref. [39]). Note that the peaks are also 
shifted slightly to the right as we deviate more from ACDM. This is also due to matter-radiation equality occurring 
at earlier times, such that only very small scales enter the sound horizon during radiation domination, therefore the 
power spectrum turnover occurs on smaller scales, also shifting the CMB peaks to smaller scales i.e. higher I. 

V. CONCLUSION 

In this paper we have investigated the effects of coupling a variable fraction of the dark matter to a scalar field 
without the presence of a potential. We find that the background evolution can deviate significantly from that of 
ACDM, depending on the strength of the coupling and the proportion of CDM coupled. This means that when 
carrying out calculations that depend upon the background, such as perturbation evolution, it is necessary to take the 
evolution of the scalar field background into account. This is due to the fact that we do not include a bare potential 
term, and so the effective potential does not have a local minimum at finite values of the scalar. For this reason, we 
also find that the late time attractor solutions often seen in quintessence models are not present here. In fact, the 
only stable solution to this model is A domination, as one might have anticipated from the cosmic no hair theorem. 

The main purpose of our work, however, was to derive the perturbation equations and solve them in order to 
calculate the CMB angular power spectrum. Here we find that there is only a small deviation from the ACDM 
predictions, even for O(0. 1) scalar coupling to matter, in gravitational units, with over half the total dark matter 
coupled directly to the scalar field. This is quite compelling and certainly a surprising result. It is, in part, down 
to the fact that the coupling itself seems to play only a small role in the anisotropics. It is actually the relative 
proportions of matter and radiation present, and the time of matter-radiation equality, that has the biggest effect 
on the CMB. The key point is that the CMB is only an indirect probe of long range dark matter interactions. CMB 
photons do not couple directly to the scalar field, and only feel its effect through their impact on structure formation. 
The scalar mediates an additional attractive force on the relevant dark matter species enhancing the growth of density 
perturbations, the effect being more pronounced the greater the coupling and the greater the fraction of coupled dark 
matter. 

It would be interesting to study the effect on alternative probes of the dark matter from supernovae to large scale 
structure, and to jointly analyse these data sets along with CMB data. As we have already mentioned, supernovae 
might be considered as a more direct probe of the long range force as it causes galaxies, dominated by dark matter, 
to follow non-geodesic trajectories. The additional attraction mediated by the scalar is also expected to increase the 
number of large clusters, which is not favoured by observation, so a detailed studied of structure formation in this 
model may also place a strong bound on our parameters. For related work with coupled quintessence scenarios and 
chameleons see e.g. Refs. [40] and [41]. 

Finally, we conclude our discussion with a word about additional long range repulsive forces. Typically these require 
a vector field to mediate the interaction, but we might have hoped to mimic its effect by reversing the sign of the 
kinetic term in our scalar field model, bearing in mind that such a toy model should not be taken too seriously because 
of violations of unitarity. However, it turns out that this trick doesn't work. The ghost-like nature of the scalar turns 
the effective potential into an unstable one, leading to rapid evolution of the scalar and time varying gravitational 
potentials. This has a dominant effect on the CMB angular power spectrum at low multipoles through a large late 
time Sachs- Wolfe effect that would not necessarily be expected in a consistent scenario involving a vector mediator. 
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